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We investigate two-species Bose gases in traps with various interactions using variational Quantum Monte 
Carlo (VMC) techniques at zero temperature. The bosons are represented by hard spheres (HS) whose 
diameter is equivalent to the s-wave scattering length in the low-energy and long-wavelength approximation. 
We explore the role of repulsive and attractive inter- or intraspecics interactions on the condensate properties 
of the mixtures, particularly the condensate fraction of each species as compared to the case when each species 
is in a separate trap of its own. We model the repulsive interactions by a hard core (HC) potential and the 
attractive interactions by a shallow model potential. The VMC density profiles and energies are evaluated 
at various interactions and two mass ratios of the species. 
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I. INTRODUCTION 

The interest in the investigation of two-species Bose- 
Einstein condensates (2BECs) in traps has grown sub- 
stantially since their first experimental realization in 1997 
[H . Since then, the literature on 2BECs in traps has ex- 
ploded II M M, 0, 0, S i, M M, M M M M, M 

[itI . [isl . fig. |2C| and the interests are now drifting more 
strongly towards these mixtures [111, , which are the 
major theme of this paper. 

In this paper we investigate the role of inter- and in- 
traspecics interactions on the properties of two-species 
Bose gases (2BEC) in a tight isotropic harmonic trap 
at zero temperature using variational Quantum Monte 
Carlo (VMC) methods. A tight trap enables us to simul- 
taneously use a low number of particles and achieve high 
densities since the volume of the trap ped cloud is much 
smaller than the usual size psi. [2^. [25| where aho ^ 10^ A. 
Here a^o = \/ h/mojho is the trap length where m is 
the mass of the boson, to the trapping frequency, and 
h Planck's constant. We represent the bosons by hard 
spheres (HS) whose hard core (HC) diameter is equiva- 
lent to the s-wave scattering length in the low-energy and 
long wave length approximation. In order to describe the 
interactions, we use a HS potential for repulsive, and a 
shallow two-body model potential for attractive interac- 
tions. Basing on this, we emphasize the qualitative na- 
ture of the present work and deemphasize comparisons 
with current experiments. The results stand alone as 
qualitative properties of the model system. A key point 
in our present research is that we do not use the scat- 
tering length in describing attractive interactions as is 
usually done in mean- field investigations, but rather the 
depth of a two-body model potential as justified later 
on. We thus vary the depth of the model potential and 
HC diameters of the bosons and investigate the resulting 
properties such as the VMC energies and density distri- 
butions. We particularly focus on the VMC condensate 
fractions and condensate density profiles in an -and to the 



best of our knowledge- unprecedented manner in the lit- 
erature concerning mixtures of Bose gases. Another key 
point here is that we focus on the factors that enhance 
the condensate depletion of the 2BEC components. We 
find chiefiy that the mixing of two Bose gases in a trap 
enhances the depletion of the condensates of each gas 
as compared to the case when either one is in a sepa- 
rate trap of its own. Thus the one-component Bose gas 
(IBEC) in our paper works chiefly as a reference system 
to which we compare our mixtures. Further we find that 
no phase separation can occur in the case of attractive 
interspecies interactions and that two Bose gases can not 
be mixed in the case of large repulsive interspecies in- 
teractions. Some of our findings are similar to those of 
Kim and Lee [2^ and Shchesnovich et al. [1] . Our work is 
particularly related to the work of Ma and Pang [13] who 
did an investigation similar to ours except that they used 
repulsive interactions only, whereas we additionally use 
attractive interactions. We further evaluate the energies 
of the systems and check them against an approximate 
model calculation. 

On the theoretical side, there have been many theo- 
ries and investigations. For example Kim and Lee [2^ 
examined the stability properties of the ground state of 
2BECs as a function of interspecies interactions. One 
of the ground states that they found had a component 
localized at the center of the trap surrounded by the 
other component thereby forming a core and shell. Ho 
and Shenoy [^ discussed binary mixtures of alkali con- 
densates and found that the heavier of the two compo- 
nents always enters into the interior of the trap and the 
lighter component is usually pushed towards the edges of 
the trap. Chui et al. [^ investigated the nonequilibrium 
spacial phase segregation process of a mixture of alkali 
BECs. Pu and Bigelow [3 presented theoretical studies 
of a 2BEC. They showed that a mixed Bose gas displays 
novel behaviour not found in a pure condensate and that 
the structure of the density profiles is very much infiu- 
enced by the interactions. 
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On the experimental side, there also have been many 
investigations. For example Modugno et al. [27j reported 
the realization of a mixture of BECs of two different 
atomic species using potassium and rubidium by means 
of sympathetic cooling. Again Modugno et al. (2^ re- 
ported on the Bose-Einstein condensation of potassium 
atoms achieved by sympathetically cooling the potassium 
gas with evaporatively cooled rubidium. Mudrich et al. 
[29| explored the thermodynamics in a mixture of two 
different ultracold Bosc gases. They showed that a hot 
gas can be cooled to a lower temperature by mixing it 
with another colder gas. Maddaloni et al. |30| demon- 
strated an experimental method for a sensitive and pre- 
cise investigation of the interaction between two conden- 
sates. They studied the effects of interaction by studying 
two completely overlapping condensates and found that 
the center-of-mass oscillations of the two condensates are 
damped if they are interacting and otherwise if they are 
noninteracting. Matthews et al. [3l| presented the exper- 
imental realization and imaging of a vortex in a two com- 
ponent BEG. They induced the vortices by a transition 
between two spin states by hyperfine splitting of *^Rb 
usinga two photon microwave pulse. Again Matthews et 
al. j32l | explored the dynamical response of a BEG due 
to a sudden change in the interaction strength and pre- 
sented a method for the creation of condensate mixtures 
using radio frequency and microwave fields. Further they 
observed an oscillatory behaviour of the condensate sizes 
when the interactions are changed. 

Although the above revealed some of the most im- 
portant properties of ultracold mixed atomic systems in 
traps, an investigation of the condensate properties and 
energies is still missing. For example, what is the role 
of the hard core (HG) of the atoms in one component in 
determining the condensate fraction of the other compo- 
nent? In previous publications [H,!!!], it has been shown 
that the hard core (HG) of the bosons plays a fundamen- 
tal role in depleting the condensate of a one-species Bose 
gas in a trap. Another issue which has not been ad- 
dressed before is the role of the mass-ratio of the bosons 
in a mixture in determining the condensate fractions and 
we briefly address this issue in this paper. 

The paper is organized as follows. In Sec. II we present 
the method we used. In Sec. Ill we outline our results 
and in Sec. IV we discuss them and connect to the pre- 
vious literature. In Sec. V we list our conclusions and 
in Appendix A we present a model for the estimation of 
the energies. 



II. METHOD 

We consider ultracold two-species Bose gases (2BEG) 
of A'^i and N2 particles, masses mi and m2, and HG di- 
ameters Oc and be, respectively, confined in a spherically 
symmetric tight harmonic trap. The total number of 
particles = iVi -I- iV2 is kept fixed and we use small 
numbers of particles since larger ones increase the com- 



putational times substantially. We investigate the 2BEGs 
using variational quantum Monte Garlo (VMG) methods 
at zero temperature. The program for VMG used in ear- 
lier publications [s^ [s^] and for the one-body-density- 
matrix (OBDM) [33] has been modified to accomodate 
2BEGs. We shall not explain the VMG technique as it 
can be found in a large number of references, rather we 
present our trial wave function and mention briefly how 
the particles are moved and how the densities are calcu- 
lated. 



Hamiltonian 



The Hamiltonian of a two-component Bose gas is 
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where mi and m2 are the individual masses of the atoms, 
Tcti, ra-iv„ are the particle position-vectors from the 
center of the trap of components a = I and 2, uji and 
UJ2 are the trapping frequencies, V\™* and 1^™* are the 
intraspecies interactions of species 1 and 2, respectively, 
and V12 * is the interspecies interaction. 



B. Units 

We take length and energy in units of the trap Cho = 
y^h/mojho and huj^o, respectively, where m = mi and 
LUho = uJi are the mass and trapping frequency of compo- 
nent 1, respectively. Using these units {H H/{HuJho) = 
H , r r/aho = r), the Hamiltonian ([1]) can be rewritten 
in the form 
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thus introducing two ratios {mi/m2) and {m2'jJ2 / miuj\) 
into the Hamiltonian. 
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FIG. 1; HCSW interatomic potential with ac = 0.05, attractive 
well range R — = 0.54 and a potential depth Vq = —3. For 
r < ac, V{r) is infinite. All lengths and energies are in trap units, 
ajio = ^/h/majfio and tu^ho, respectively. 



C. HCSW interactions 

We model the two-body interactions by using a hard- 
core square weU (HCSW) potential. Essentially, it is a 
hard core plus an attractive tail added to it. Figure [T] 
shows our model potential where V{r) is the depth and 
r the two-body intcrparticlc distance, all in units of the 
trap. Here for example the bosonic HC diameter is Uc = 
0.05, the depth is Vb = —3, and the range is d = R ~ 
Uc which we keep fixed at 0.54. This range is of the 
same order of magnitude as that used in a previous work 



[35l | for another potential and we return to this point 
in Sec. IIVDI In this paper we are chiefly interested in 
using the depth of the HCSW to describe the interactions 
and not the associated scattering length. Nevertheless we 
check the stability of the systems at the first Feshbach 
resonance when a — s- zLoo in Sec. IIV El later on. 



D. Trial wave function 

The general form of the trial wave function is 

*T({ri},{ai};{r2},U22};te2}) = 

n n /i2(ii"i™ - ""ini) X 
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where {r^r} = {vai, ...^r^N^), goijcri) are single particle 
wave functions for particles of type a = 1,2, and fa^a^ 
are pair correlation functions for intraspecies and inter- 
species interactions with variational parameters given by 
the sets {^critra} = {/3<ti<t2 , Tcriaa , etricrs }■ Hcrc the pairs 
{ci, CF2} are {11} and {22} for intra- and {12} or {21} for 
interspecies interactions. There can be several choices for 
the Jastrow functions depending on the interatomic in- 
teractions. In our case we constructed a flexible Jastrow 
function inferred from the exact solution of two particles 
interacting via a HCSW potential, 
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where ro-io-aO is the position of the maximum of the 
Jastrow function and aaia2 is the HC diameter, where 
ail = flc, 022 = be, and ai2 = a2i = {ac + hc)/2. The 
sinusoidal part of ^ is taken similar to the exact solu- 
tion of two particles colliding inside a HCSW with rel- 
ative energy Ea^„,_ = ^^Kia^li 2^0-1 0-2) and HC diame- 
ter aaicT2 by replacing the HCSW wavevector K 0-^0-2 = 
•\/2/Xo.i(T2 (Kricr2 + £'cri(T2)/^^ for each type of interaction 
of strength Va--^a-2 by a variational parameter 70-1 0-2 ■ This 
is in order to decouple the Jastrow functions from their 
HCSWs and to introduce some flexibility to them. An- 
other reason for this replacement is that we do not know 
the values of E^-^a-2 at the higher densities and we there- 



fore allow 7o-iCT2 to vary slightly in order to indirectly 
imply a value for i?cri(T2- (Here /in = mi, /i22 = 
and fii2 = mim2/{mi + 1712).) In our simulations, 
the optimized 7^10-2 is always very close to A'cri(T2 and 
7criCT2 ^ Ka--^a2- Thus attractive interaction between 
the particles is included in the Jastrow via 7,71 0-2 ■ We 
then join the sinusoidal part at Vij — r^-^a-20 to another 
function which decays to 1 in the long range. Note then 
that ro-io-ao is not necessarily equal to Rai<72: the edge 
of the HCSW for each interaction type, and depending 
on the well depth it can be either inside or outside the 
HCSW. The reason for this construction is to provide 
smooth Jastrow functions whose maxima are at inter- 
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particle distances large enough to bring the bosons close 
together. Further it is important to note that the attrac- 
tion between the bosons is mainly caused by the Jastrow 
function ([4]), particularly by the "bump" of the Jastrow 
which is higher than 1 at r = ra^cr^ii- The part of the Jas- 
trow function in the range aaia2 < fij ^ t 01020 is then 
repulsive, whereas in the range > r„-^„^Q attractive. 
Note that in the case of only repulsive interactions the 
HCSW depth Va^a^ and /S^icra are set to zero. Therefore 
when 7^^^2 0, 

7cricr2~^0 ■7(71(72 ^ij \ ^ ij / 

brings us back to the HS Jastrow function. Ao-io-a (and 
''cri(T2o) E^re parameters that join the Jastrow in the 
interparticle-separation range a„^u^ < < r„-^„^Q to 
that in the range ry > r„-^cr-20 foi' the same slope and 
amplitude. For the single particle wavcfunctions, wc use 
Gaussians of the form 

g„{v„i) = exp(-Q;<jr^J (6) 

where are variational parameters. Later on in this 
paper we shall see that even with Gaussians centered at 
the origin, the variational wavefunction ([3]) is still able 
to describe phase separation. This indicates that the 
trial wavefunction is dominated by the pair correlation 
functions /o-icrj rather than the single particle functions. 
In fact the Gaussians chiefly cause the density to vanish 
at the edges of the cloud thus indirectly confining the 
cloud within a certain volume. 



E. Moving the particles 

The particles are moved according to 

= r,, + Ar, • (77 - 0.5) (7) 

where v'^^ are new positions, 77 is a random number be- 
tween and 1 and Arc- are step size vectors which are ad- 
justed to obtain optimal diffusion through configuration 
space- i.e. to obtain a VMC acceptance rate of ~ 50%. 
After each update of the particle coordinates the pro- 
posed move is either accepted or rejected according to 
the (MRT)^ [3^ algorithm where the square of the trial 
wave- function is used as the probability distribution from 
which particle configurations are sampled. 



The variance is given by a\ = {E"^) — {E)"^ where in 
general for any operator O 

with vl/o = *({ri},{C0i};{r2},U02};U?2}) and 
vPi = *({ri},{eii};{r2},U22};Ui2}) and /dn = 

nf=\/rf'^i., Sdv, ^ Y{%Jd^r,,. Here |vfi|Vl*oP 
are weights used for the reweighting process of the vari- 
able O. {^^1(72} ^'^^ initial and {(,0102} the optimized 
sets of variational parameters. In the Gaussians we also 
use the initial q;° and the optimized a^- 

G. Condensate fraction 

In the systems considered here we have two condensate 
fractions n"^^ and rS^^ for components 1 and 2, respec- 
tively. The overall condensate fraction of the mixture is 
{n^^Nx + n^o^^N2)/{Ni + N2) but we only focus on the 
individual rig^'' and tt.q^'' . The condensate fraction of each 
component is evaluated by calculating the eigenvalues of 
the natural orbitals using the one-body-density-matrix 
(OBDM) of each component in a manner similar to a 
calculation by DuBois and Glyde js^l . By using the trial 
wave function of Sec. Ill Dl wc evaluate the OBDM for 
components 1 and 2, respectively, as follows. In order to 
make the equations more compact, we define 

'^12, ■ • ■ , riN, , {61}; {i"2} , {£,22}; {62}), 

ipir[j^,ri2, ■■■ , rijvi , {61}; {^2} , {£22}; {62}), 
Q2{r2i,---) = 

, {Cii}; ''21, '"22, • ■ • , ^2^2, {62}; {62}), 

Q2{r'2,r--) = 

, {Cll}; ?'21) ''22, • • • , r2N2,{£22}; {^.u})- 

Hence the OBDMs are written 



F. Minimization of the variance of the energy n (r r' \ — 

In order to optimize the trial wavefunction, we numeri- Y[i=2 I ^ 'riiY[j=i J d r2jQi{rii, )Qi(''ii7 ) 
cally minimize the variance of the energy cte with respect < ip\^p > 

to the variational parameters of the trial wavefunction. (9) 
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leaving out the integration over rn and and 



A. Stability of the mixtures as compared to 
one-species Bose gases 



Pb{r2i,r'2i) 



< I V' > 



leaving out r2i and ■ Here 



(10) 



<^\^>= J dn y dr2|V;({ri},Ui};{r2},{6};Ui2})P, 

is the normalization factor and Citi<72 ^^'^ the optimized 
variational parameters. Hence we extract the OBDM 
for each component from the two-body density matrix 
(TBDM) of the mixture by integrating out the contri- 
bution from the other component. In a manner simi- 
lar to Ma and Pang [l^ then, each component is es- 
sentially treated as a subspecies with its own properties 
but still it is not independent of the other species as a 
result of the interspecies interactions. The interspecies 
interactions arc included in the OBDM through the in- 
terspecies Jastrow function /12. From the trial wave func- 
tion wc can verify that Eqs. ([9]) and (fTO|) reduce to the 
one-component case if the interspecies interactions are 
turned off. That is the interspecies Jastrow function /12 
becomes equal to 1 and the two components become in- 
dependent of each other as they are now noninteracting. 



H. Density profiles 

The densities are calculated during a VMC run by di- 
viding the space along the radial direction into spherical 
shells (bins) concentrated at the center of the trap and 
collecting the particles of each species in them as was 
done before by DuBois and Clyde [s^. 



III. RESULTS 

In what follows we present the results of our Monte 
Carlo simulations. We display and discuss the resulting 
VMC density profiles and the condensate fractions of our 
mixtures with various interactions. We further compare 
the condensate fractions of the mixtures with the con- 
densate fractions of their components when either one is 
in a separate trap of its own. We compare our VMC en- 
ergies with the results from an approximate mean-field 
model derived in Appendix A. Wc further reveal the role 
of the mass ratio toi/to2 hi determining some properties 
of the Bosc gases. The trapping frequency is set to be 
the same for both components (wi = U2) and the mass 
ratio is arbitrarily chosen to be mi/m2 = 1.200. 



During the numerical optimization of the variational 
parameters as explained in Sec. IITfI we plot the energy 
EvMc/N versus the set of variational parameters {£,ai<T2 } 
used in our wave function. The numerical optimization 
procees changes the variational parameters over several 
iterations and searches for a minimum in the energy vari- 
ance which also leads to a minimum in the energy. After a 
number of iterations, we obtain plots such as those shown 
in the following figures. For example Fig. [5] displays the 
VMC energy Evmc /N vs one of the variational parame- 
ters ai for a mixture with ^ O.l, ~ 0.2 and repulsive 
HC interactions only (upper frame) and for a mixture 
with Qc = 0.2, be = 0.3, attractive (HCSW, V12 = -10.0) 
inter- and repulsive HC intraspecies interactions (lower 
frame) . The figure depicts clearly the presence of energy 
minima at ai ~ 0.25 and ^ 2.9, respectively. The behav- 
ior of the energy versus the other variational parameters 
is the same as in Fig. [2] and all of them display energy 
minima. After the completion of the optimization pro- 
cess and in the final evaluation of the wavefunction for 
each system, wc choose the variational parameters that 
correspond to the energy minimum, i.e., the ground state. 
All of our repulsive or attractive 2BECs display energy 
minima as above and we can therefore state safely that 
our mixtures are stable systems. 

In comparison, Fig.[3]shows the VMC energy against a 
for a HCSW IBEC of 20 particles, HCSW depth V = -6 
and Oc = 0.2 using the same trial wavefunction ^ but 
set for one component only. The figure shows a peculiar 
result, namely the presence of two equal energy minima 
at a ~ 1.15 and 1.2, i.e. a degeneracy. One of the minima 
is due to the single-particle, the other due to the Jastrow 
part of the trial wavefunction. The single-particle wave- 
function is connected to the external trapping potential 
and the Jastrow function to the interparticle interactions 
and generates as such the energy minima due to these 
potentials. This plot has been generated from two VMC 
runs using different minimization directions 37[ in order 
to ensure the presence of the two minima. We have seen 
this phenomenon in all the VMC runs for this particular 
system at various other HCSW depths. We do not un- 
derstand at the present why this double minimum does 
not occur in 2BECs. 



Definitions of the densities 



In order to describe the density of the systems, we used 
na^ with Oc the HS diameter of the bosons for compo- 
nent 1 and with be the HS diameter for the bosons 
of component 2. We thus describe the systems by the 
HC density only, even in the presence of attractive inter- 
actions. We define the total VMC spacial density distri- 
butions (condensate + normal parts) by 711(7") and 712(7") 
for components 1 and 2, respectively, and in units of af^^ 
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display the properties of both components as a function 
of their HC densities in a single plot. For this particu- 
lar purpose, we use a unified term, namely na^g with 
qhs the HS diameter Oc or 6c, to describe the HC den- 
sity of either component at the center of the trap and 
na^g is used then under the following conditions. In a 
single HC or HCSW IBEC 7ia%g = n{0)al where 7i(0) 
is the number density at the center of the trap and Uc 
is the HC diameter of the single-species bosons. Since 
we are dealing with more than one species, na^c; of each 
component has to be defined for various cases of inter- 
actions in the 2BECs. For systems with attractive inter- 
species interactions where there is full mixing (Sec. IHI C|) 
na^g ~ ?ii(0)a^ for component 1 and na^g = n2(0)6^ 
for component 2. In the case of repulsive interspecies in- 
teractions na'fjg is the HC density of the core only since 
the two species phase-separate (Sec. IIIIDI and lIIIE[) and 
it is difficult to define na^g for a shell. In all of our in- 
terpretations, we do not consider na^g to be the overall 
total density (nxir) — ni{r) + n2(r)) of the mixture at 
any time. 



Attractive inter- and repulsive intraspecies 
interactions 

1. Density profiles 



FIG. 2: VMC energy per particle versus the variational parameter 
ai of the trial wave function (|3]l for a trapped Bose gas mixture 
of A^i = 20 and N2 = 10 particles. Upper frame: HC 2BEC with 
ac = 0.1 and be = 0.2, lower frame: HCSW 2BEC with = 0.2, 
be = 0.3 and V12 = -10.0. V12 is the interspecies HCSW depth. 
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FIG. 3: As in Fig. [2]but for a HCSW IBEC with Af = 20 particles. 
The depth of the HCSW is V = —6 and the arrows indicate the 
locations of two equal energy minima {E/N = —1.740). This plot 
is a result of two VMC runs of the same system. 



where r is the distance of a boson of either species from 
the center of the trap (r = 0). Correspondingly ri,o,i(r) 
and no,2{r) arc the VMC condensate density distribu- 
tions. The total VMC density of a IBEC is written n{r). 
In the interpretation of our results, we sometimes need to 



The goal of this and the following sections is to display 
the spacial VMC density profiles of 2BECs with various 
interaction parameters which have largely not been dis- 
played before in the literature. Figure |4] displays the 
MC density profiles of 2BECs with attractive (HCSW) 
interspecies and repulsive (HC) intraspecies interactions. 
The points with error bars represent nx^r) of the mix- 
ture, the open circles represent the density 711(7") and the 
open triangles 712(7"). The thick and thin dashed lines 
are 779,1 (7^) and 770,2(7"), respectively. The strength of the 
interspecies interactions is indicated by the depth of the 
HCSW, \Vi2\, and the range of the HCSW is kept fixed 
at 0.54. In all our mixtures, here and thereon, the com- 
ponents have Ni = 20 and N2 ~ 10 particles and in this 
section ac = 0.2 and be ~ 0.3. There is a particular rea- 
son for the choice of the latter large Uc and be above, this 
is in order to enable substantial depletion of the conden- 
sates. We keep ac and be fixed and increase V12 from 
to —40 in the "negative" sense. 

The key features of Fig. 5] are as follows. The attrac- 
tive forces enable full mixing of the two components. In 
frame (a) component 2 is slightly pushed out towards the 
edges of the trap due to the repulsive interspecies inter- 
actions arising by setting V12 = 0. That is when the 
attractive part of the HCSW is switched off, the HCSW 
changes to a repulsive HC potential. Then we note that 
although V12 = 0, full mixing of the two components 
is still possible. At the instant the HCSW is "switched 
on" as in frame (b) , component 2 is pulled back towards 
the center of the trap with no remnant expulsion at the 
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FIG. 4: VMC density profiles and condensate properties of 2BECs with HCSW inter- and HC intraspecies interactions. V12 is the 
interspecies HCSW depth. Points with error bars: nT{r), open circles and triangles: ni(r) and n2(r), respectively, thick and thin dashed 
lines: corresponding nQ^i(r) and riQ 2(''), respectively. 



edges of the trap. The densities in frame (b) jump now 
above those in frame (a) and continue to rise as V12 is 
increased. In frames (a) and (b) the condensate densities 
no,i(r) and no,2('') are similar in shape to their corre- 
sponding total densities ni(r) and n2{r), but in the rest 
of the frames (c-f ) they are not. Rather they obtain a flat 
shape in frames (c) and (d) after which they are slightly 
pushed out towards the edges of the trap in frames (e) 
and (f). In any case, the attractive forces prevent the 
condensate from total expulsion towards the edges of the 
trap. We anticipate that as the density rises further with 
V12, the condensates wifl be pushed out further towards 
the surface of the cloud because the condensate seeks the 
lower density regimes of the cloud. The reason is because 
the lower cloud density at the edges of the trap causes a 
lesser local condensate depletion than the higher density 
towards the center. Note also that the total densities in 
frames (a-c) have a Gaussian shape, but then they divert 



from it somewhat. The densities rise also significantly 
with the increase of attractive interspecies interactions: 
nT{r) rises by a factor of w 18 from frame (a) to frame 
(f) as the cloud radius shrinks in size by a factor of w 3. 

We are also able to use large attractive interspecies 
interactions {V12 = —40) at the energy scale of ultracold 
Bose gases and still obtain energetically stable systems. 



2. Condensate fraetions 

The goal of this section is to display the effect of com- 
plete mixing on the depletion of each condensate in a 
2BEC mixture as compared to the case when cither con- 
densate is in a separate trap of its own in which case it 
forms a IBEC. For this purpose we consider the mixtures 
in Fig. m Figure [5] displays their condensate fractions as 
a function of na^g. The open and solid squares represent 
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and solid triangles: components 1 and 2 of the HCSW 2BEC and 
Vi2 is the depth of the HCSW for some of the points. The points 
are larger than the error bars. 



the condensate fractions no for HC IBECs with N — 20 
and = 10 particles, respectively, which act as our ref- 
erences. The open and solid triangles display the conden- 
sate fractions of the mixtures n^"^ and n.Q^\ The figure 
depicts clearly that the depletion of the condensates in 
the 2BECs is larger than the IBECs. This reveals that 
mixing enhances the depletions of the constituent con- 
densates due to their interspecies interactions. A signif- 
icant feature is that tIq^^ and n^^ as a function of na^g 
coincide at the larger nuj^g. We may attribute this to 
the fact that since the two components are completely 
mixed the system behaves similarly to a IBEC. 



D. Attractive intra- and repulsive interspecies 
interactions 



1. Density profiles 

Figure [S] displays density profiles as in Fig. 2] but with 
Gc = 0.2 and be = 0.4 with HCSW intra- and HC inter- 
species interactions (V12 = 0). We keep Uc and be fixed 
and vary the intraspecies HCSW depths {Vu, [i = 1,2)) 
in the range Vu = —4 to —16 keeping Vu = V22. As a 
result of the attractive intraspecies interactions, the core 
and the shell contract in volume and the density of the 
system grows substantially in response from nT(0) ~ 2.5 
in frame (a) to ^ 14 in (d) as the radius of the cloud 
shrinks from ~ 3 to 1.5. The shell is pushed radially 
inwards towards the center of the trap by the confining 
forces of the trap. Contrary to Fig. 21 the no,i(r) and 
nQ^2{f) profiles keep following the shape of their corre- 
sponding ni(r) and n2(r). 



2. Condensate fractions 

In what follows we investigate the condensate prop- 
erties of 2BECs, this time with attractive intra- and re- 
pulsive interspecies interactions. This is somewhat the 
opposite case of Sec. IIII CI where attractive inter- and 
repulsive intraspecies interactions are used. We chiefly 
aim at revealing the difference in the results when using 
different types of combinations of repulsive and attrac- 
tive interactions. Figure [7] compares now the condensate 
fraction rig^'' of the 2BEC of Fig. [6] (solid triangles) with 
HCSW intraspecies interactions (intra.) against nQ^-* of 
the 2BEC of Fig.g] (open triangles) with HCSW inter- 
species interactions (inter.). The open circles display the 
condensate fraction uq of a HCSW IBEC of = 20 par- 
ticles and Qc = 0.2 in a separate trap of its own, this is 
our reference system. Here the HCSW depth is varied in 
the range {V ~ to —16). 

The condensate fraction n^^'' of the mixture with 
attractive interspecies interactions (inter.) shows the 
largest depletion, although the bosons in its shell have a 
smaller HC diameter than those in the mixture with at- 
tractive intraspecies interactions (intra.). It seems that 
the attractive intraspecies interactions boost the value 
of the condensate (solid triangles) beyond the HC in- 
traspecies result (open triangles). It will be shown in Sec. 
IIII El below that for a number of purely repulsive mix- 
tures with the same bosonic HC diameters in the cores 
but different HC diameters in the shells, the condensate 
depletion is larger in the 2BECs with larger bosons in 
the shell. Our reference shows again the smallest con- 
densate depletion which is again a manifestation of the 
fact that mixing enhances the condensate depletion be- 
yond the IBEC result. 

Substantial depletion is observed in FiglHlat Vu — —16 
(^ 50% for ng^'' and ~ 60% for rig^^). If compared to 
Fig. d] we can see that this amount of depletion sets in 
there at V12 = —15 (frame e), that is at a comparable 
(interspecies) HCSW depth. However, the density in the 
latter 11) is lower than the former {^-^ 14) because 
the number of attracting pairs of two-species bosons is 
smaller in Fig. 3] than in our case here which leads to a 
slower rise in the density with HCSW depth. 



E. Repulsive inter- and intraspecies interactions 

1. Density profiles 

Figure [8] displays the VMC spacial density distribu- 
tions of HC 2BECs with various interspecies interactions 
and the same definitions of points as in Fig. U) Here Oc 
is varied while be is kept fixed at 0.3. We can see that 
mixing of the two components is enabled up to Oc = 1.0 
before complete phase-separation sets in. On increasing 
flc beyond be in frames (a-h), component 1 (of HC Oc) is 
gradually pushed out towards the edges of the trap with 
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FIG. 7: Condensate fraction tIq^^ for HCSW 2BECs compared to 
a reference. Open circles: (reference) HCSW IBEC witli TV = 20 
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tween 2 and 16. Solid triangles: 2BEC of Fig. |6] with attractive 
intraspecies interactions. Open triangles: 2BEC of Fig. |4] with at- 
tractive interspecies interactions. The numbers between brackets 
near to some of the points show \ Vij\ for the corresponding systems . 
The points are larger than the error bars. 



and LUi ~ LU2. This can be explained as follows. Essen- 
tially, as the HC diameter of the bosons of one component 
increases, the Bose gas expands in size in order to acco- 
modate the larger bosons. As a result, the component 
with lower intraspecies repulsion "falls" into the center 
of the trap seeking the minimizaton of the total repulsive 
potential energy. As the HC diameter Oc is increased, the 
interspecies repulsion rises pressurizing the core radially 
towards the center of the trap. As a result, the density 
of the core n2 (0) rises by a factor of « 2 from frame (a) 
to (h) and nearly as from frame (d) on, ni(0) begins to 
approach zero and the two components begin to sepa- 
rate into a shell and a core. If we imagined removing the 
shell completely from the trap, the core will expand and 
become almost uniform in density and therefore "flat" in 
shape [H, [U • We found that it is very hard to "squeeze" 
the core further to higher density by increasing Uc beyond 
2. In frame (h) total phase separation has occured leaving 
a dip at the boundary between the two components. One 
could imagine placing a third-species particle in that dip 
as it is a potential trap by itself. Contrary to the case of 
Figs.HJd-f), the condensate density distributions of both 
components follow the shape of their total densities up 
to phase separation. 



the rise of interspecies repulsion [qc + bc)/2. This is con- 
trary to our expectations because we thought that com- 
ponent 2 (of HC be) would be pushed out instead since 
it is the hghter of the two given that mi/m2 = 1.200 



With repulsive interactions only in these mixtures, we 
can always have stable systems if there is sufHcient re- 
pulsion between the bosons of the core counteracting the 
outer pressure arising from the shell. Otherwise a dilute 
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core collapses readily under the heavy pressure of a dense 
shell. 

Figure [5] displays the MC density distributions of a 
2BEC with HC interactions only where Oc and be are both 
increased at a fixed ratio ac : be = 1 : 2. In this case the 
density of the core ni(0) decreases because both ac and 
be are increased. A peculiar result is that even at very 
large values of the interspecies repulsions no complete 



phase separation is observed as it occurs in Fig. [8ji. Some 
uniformity in the density distribution of the core arises 
at the larger ac- 

In Fig. [TU] wc make comparisons between densities at 
the center of the trap as a function of the interaction 
parameter NiUc for various systems with repulsive inter- 
actions. The crosses are for a HC IBEC of 20 particles 
in a separate trap of its own and the same trap length as 
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FIG. 9: As in Fig. |8]but with the HC diameters both varying at a fixed ratio of ac : be = 1 : 2. 



before a^o = y^h/niiuji. The open circles and triangles 
are, respectively, for the core in Fig. [9] and an additional 
mixture with Oc : 6c = 1 : 3 whose density profiles we do 
not reveal. The solid diamonds are for the core in Fig. [S] 
Thus the goal is to show the effect of mixing a HC IBEC 
with various other IBECs on the central core density of 
the system. The density n{0) of the HC IBEC drops as it 
is mixed with another component, and for a larger ratio 



of be relative to Oc the core density drops further. The 
density of the core with one of the HCs fixed {be = 0.3) 
varies only slightly as Qc increases. Note that the values 
of the core density at r = 0, except for the latter case, 
converge at the higher NiUc- 
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FIG. 10; Density at the center of the trap (r = 0) versus the 
interaction parameter N\ac for HC systems only. Crosses: 
IBEC with Ai" = 20 particles {reference), open circles: 2BEC 
of Fig. [9] with Oc : be = 1 : 2, open triangles: 2BEC with 
flc : 6c = 1 : 3, solid diamonds: 2BEC with be fixed and 
varying (Fig. 15)). The points are larger than the error bars. 



FIG. 12: Condensate fraction versus na^g for HC systems only. 
Open squares: IBEC of A'^ = 20 particles, open triangles and cir- 



cles: riQ^' and tiq'' for the 2BEC of Fig. (9] solid triangles and 
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FIG. 11: Condensate fraction versus na^g for a HC IBEC and 
component 2 in 2BECs. Open triangles: HCSW 2BEC of Fig.E] 
solid triangles: HC 2BEC of Fig.|8] open circles: HC IBEC with 
N = 10 particles (reference). The points are larger than the error 
bars. 



2. Condensate fraction 



Figure [Tl] compares the condensate fraction Hq of- 
Fig.[5] (oc : 6c = 1 : 2) to two other systems. The open 
circles display no for a HC IBEC of 10 particles in a 
separate trap of its own and the same trap length as be- 
fore Oho ~ y^h/miLUi. The solid triangles display Uq^' 
for the HC 2BEC of Fig. [5] with be fixed and Oc varying. 

The open triangles display n*,^^ for the HCSW 2BEC of 
FiglHwith Qc and be fixed and V12 varying. We note that 
the condensate depletion is highest in the HC 2BEC of 
Fig. [SI The depletion of the condensate in the HCSW 



2BEC is less pronounced. This reveals that the repulsive 
interspecies interactions play a more pronounced role in 
depleting the condensates of the mixture than attrac- 
tive interspecies interactions. The attractive interspecies 
interactions boost the condensate somewhat above the 
HC-interspecies interactions result. The depletion of the 
condensate is lowest for a IBEC. Thus the mixing of con- 
densates enhances their depiction due to the presence of 
interspecies interactions of various strengths as compared 
to the case when they are separate, each in a trap of its 
own. 

Figure [12] displays chiefly the condensate fractions of 
the components as a function of na^ at the center of 
the trap in various HC systems compared to a reference. 
The open squares represent the condensate fraction uq of 
a HC IBEC of 20 particles {reference). The open trian- 
gles represent the condensate fraction n'^^ of the 2BEC 
of Fig. [9| where Oc and be are varied at a constant ratio 
of 1:2 and the solid triangles that at a ratio of 1:3, re- 
spectively. In addition, and for further comparison, the 
open and solid circles represent Uq of the latter two 
mixtures, respectively. We can see that n^^^ is lower for 



1 : 3 than for 1 : 2 and the same is true for 



.(2) 



The latter values of Uq are lower than those of the HC 
IBEC displayed for comparison. This shows again that 
mixing and a larger interspecies interaction enhance the 
depletion of the condensates in each component beyond 
the IBEC result. 



F. Energies 

In this section we compare our VMC energies for HC 
2BECs against the energies calculated by an approximate 
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model derived from mean-field results in Appendix A. 
The estimate that we obtained for the total energy of 
the mixture E'^'^^ is given by Eq. (jAlOp where MFA 
stands for mean-field approximation. 

Figure [13] displays our VMC energies for the repulsive 
mixtures investigated in Figs. [9] with Qc '■ be = 1 : 2 (open 
triangles) and [8] with be = 0.3 (solid circles). The open 
diamonds show the additional VMC calculations for HC 
2BECs with Oc : 6c = 1 : 3 similar to Fig. [9l The crosses, 
times, and open circles display E^'^^ for the systems 
indicated by the open triangles, diamonds, and solid cir- 
cles, respectively. We note that there is good agreement 
between the energies £;*^/^ and the VMC results at the 
lower but then they begin to diverge somewhat at 
the higher densities. The same is true for iJ^g^"^ . E^'^^^ 
largely do not agree with the VMC results but show the 
same trend in their values. This might be chiefly due 
to the fact that the TF radius of the core is not a good 
representation of the cloud radius for this particular case 
because the two components completely phase-separate 
at the higher interspecies repulsions. In all cases, the 
energies rise with the HC densities. The rise is steepest 
when he is fixed and Oc varied, the reason being due to the 
fact that the core density varies slowly with the rise of Cc 
(see Fig. IS]). That is as the shell is expelled towards the 
edges of the trap, the potential energy rises faster than 
the change in ^2(0)6^ compared to the other systems. 



G. Effect of mass ratio 

In this section we present the role of the mass ratio 
"niratio = mi 111112 m determining the properties of mixed 
Bose gases. We consider two mass-ratios, the previous 
"rriratio = 1-2 and a new rriratio = 5 and in order to keep 



the trap length unchanged, we only change m2. We fur- 
ther note that changing rriratio changes the energy of the 
system since mi/m2 and m^lrtix appear explicitly in the 
Hamiltonian ([2]). 

We consider the systems depicted in Fig. [9] with 
TTiratio = 1-2 and Compare its properties with those of 
exactly these same systems evaluated at rriratio = 5.0. 
Figure [14] displays the condensate fractions of the latter 

systems vs na^ where the condensate fraction rig^^ (and 
(2) 

Uq ) is the same for both values of rriratio- Therefore 
rriratio has no influence on the relation between conden- 
sate fraction and HC density. The scenario is however 
different if one plots the condensate fractions as a func- 
tion of the HC interaction parameter NiUe as in Fig. 

1151 The crosses and open circles show and the open 

(2) 

and solid triangles each for rriratio = 5 and 1.2, re- 
spectively. We can see that the condensate fractions for 
rriratio = 5 are higher than for 1.2 because the central HC 
densities are lower for 5. Effectively, as m2 is reduced to 
increase rriratio, the trapping forces confining the shell 
{—V\m2to'^r\) are reduced accordingly. Thus, the cloud 
of the mixture expands as the pressure of the shell on 
the core is lifted causing the central densities to decline 
in favor of an increase of the condensate fractions. Fig- 
ure [T6| displays na^ as a function of iViOc for the latter 
systems where na^ for rriratio = 5 is lower than for 1.2 at 
higher NiUe as explained above. As we increase rriratio, 
the energy of the system rises as demonstrated in Fig. 
1171 where Ey^c/N is plotted against naj? for two rriratio 
values and Oc : 6c = 1 : 2. 
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FIG. 17: VMC energy for the same systems in Fig. 1161 and two 
mass ratios. The points are larger than the error bars. 



IV. DISCUSSION 

We discuss now important details on the results of our 
calculations and connect to the previous literature. We 
elaborate on the mixing and phase separation of com- 
ponents, stability of the mixtures and the origins of the 
enhanced condensate depletion. First of all, however, we 
mention briefly the work of Ma and Pang [l3| which is 
most relevant to ours. 



A. The work of Ma and Pang 

Ma and Pang [l3l investigated HC 2BECs trapped in 
a three-dimensional isotropic trap at finite temperature 
using path-integral quantum Monte Carlo simulations. 
Their main interest was in the structure of the mixtures, 
i.e., the densities and their profiles. They particularly 
concentrated on the conditions under which phase sepa- 
ration occurs and treated the two-species system as two 
subsystems each of which contains one species with its 
own statistics. However, the two species do not behave 
independently of each other due to the interspecies in- 
teractions. They found that by changing the mass ra- 
tio of the components 1122/ mi the lighter particles are 
pushed outward and form a shell surrounding the heav- 
ier core. Further, the density and condensate fractions 
of the mixture drop with a rise of the interspecies inter- 
actions. When identical external potentials are used, no 
phase separation is observed, but when they are different 
phase separation occurs. They also found that the spa- 
cial phase separation is independent of 7712 /mi and that 
the species with the larger scattering length favors the 
formation of a low-density outer shell. 

Now in our work here we conducted our calculations 
at zero temperature and we simulated two-species Bose 
gases with both attractive and repulsive interactions. We 
particularly concentrated on the role of the inter- and 
intraspecies interactions on the enhancement of the con- 
densate depletions in the mixture as compared to the 
case when each component is in a separate trap of its 
own. We followed Ma and Pang in treating each compo- 
nent as a subsystem with its own properties. In addition 
to their investigations on the effect of the mass ratio, we 
additionally investigated its effect on the energies and 
condensate fractions. 



B. Mixing and demixing 

In the case of intercomponent attraction as in Fig. [4] 
no phase separation occurs because the two gases attract 
each other. Note that the attractive interspecies interac- 
tions at some point overwhelm the repulsive intraspecies 
interactions as identified by the large increase in the over- 
all density nxir) of the system. That is, the repulsive 
intraspecies interactions do not become attractive, they 
just get overwhelmed similar to a case discussed by Cliui 
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and Ryzhov [38|. Further, the two gases are now trapped 
by the attractive potential of each other such that the 
importance of the external trap is underplayed. Since 
in all our calculations the minima of the two confining 
potentials coincide, the two components interpenetrate 
completely [l2j and are drawn together into what resem- 
bles a IBEC system acting similarly to it. We conclude 
that in this case, the attractive interspecies interactions 
play a more pronounced role than the repulsive ones in 
determining the properties of the systems. Therefore the 
condensates are able to migrate towards the edges of the 
trap at the higher densities na^g as in the case of a HC 
IBEC [11]. In the somewhat opposite case of Fig. [6] the 
intercomponent repulsion pushes the condensate of the 
core towards the center of the trap. 

In the case of Figs. [5] and [51 we find chiefly that full 
mixing is impossible at large repulsive interspecies inter- 
actions. Larger interspecies repulsion leads to full phase 
separation as in Fig.[5h and even though the shell there 
has a much lower density than the core, substantial de- 
pletion ~ 20% is still observed which is attributed to the 
presence of large bosons in the shell. Note that it is hard 
to define a density na^g for the shell as it is expelled 
towards the edges of the trap. 

An investigation of the detailed nature of the over- 
lap region between the shell and core is also important 
since it influences properties such as the ground state 
energy, the excitation spectrum and the coUisional relax- 
ation rates as outlined earlier by Barankov [ll|. He ex- 
plored the boundary between two repulsively interacting 
condensates in the weak and strong separation limits and 
found that the asymptotic behaviour of each condensate 
far from the boundary is determined by its correlation 
(healing) length. In the case of strong separation, he 
found that there exists a hollow in the total density pro- 
file which is very deep. The latter allows the investigation 
of one-particle excitations at the boundary between the 
components as well as surface wave excitations due to the 
surface tension [ll|. As a result of the full separation in 
Fig. [U we also observe a hollow in the total density profile 
between shell and core as discussed by Barankov. In the 
future one could add one foreign particle to be trapped 
by this hollow and investigate its energy as a function 
of some property of the mixture using the Monte Carlo 
method. Such a hollow is however not observed in Fig.[S] 
because no complete phase-separation occurs. The rea- 
son is because as the HC diameters of the bosons in the 
core are increased, the bosons spread out and the core 
expands. As a result, these bosons penetrate into the 
shell which is pushed in the opposite direction towards 
the center of the trap by the confining forces of the ex- 
ternal potential. This is however not the case in Fig. [5] as 
the HC diameters of the bosons in the core are kept fixed 
and that of the shell increased. Thus complete phase- 
separation is only possible when the size of the bosons 
in only one component is increased. This has also been 
confirmed by Ma and Pang previously. 

Two length scales can be used to characterize a two- 



component BEC [T^: one is the penetration depth, the 
other is the healing length. The penetration depth is a 
measure for the width of the overlap region and, as we 
can see from above, a function of the interspecies inter- 
actions. It is largest in the case of attraction between the 
two components as it is equal to the radius of the cloud 
whereas in the case of large intercomponent repulsion it 
is smaller than the radius of the cloud. The penetration 
depth is reduced as the intercomponent repulsion rises. 
Nevertheless, complete mixing is still possible at moder- 
ate repulsive interspecies forces as demonstrated in Figs. 
[5] and [HI Since the density profile of the core in the latter 
two systems is very much influenced by the presence of 
a shell, we anticipate that the healing length of a IBEC 
changes upon mixing with a shell. In what follows, we 
discuss some of the previous literature in connection to 
our current observations. 

Shchesnovich et al. 01 studied a ^^Rb and ^^Rb 2BEC 
by varying the interspecies interactions. They found that 
these two components would not separate if the inter- 
species interactions are attractive and in this paper we 
have verified this point as well. They argued that a sep- 
aration of the two species takes place when the energy 
gain due to the attractive intraspecies interactions over- 
whelms the quantum pressure at the interface of the two 
species. 

Cornell et al. [T^j reviewed some early results on 
mixed condensates and provided a qualitative exegesis 
of the theoretical and experimental techniques that are 
involved. They found that there is a critical value for the 
interaction term a^2 = \/aiCL2 beyond which phase sepa- 
ration occurs with little spatial overlap. This is when the 
scattering length ai of component 1 becomes larger than 
02 of component 2, causing atoms 1 to move favorably 
towards the edges of the cloud forming a spherical shell 
around the core consisting of atoms 2. Our results are 
in line with those of Cornell et al. [l2| and also Hall et 
al. [l3| as we also observe that the component with the 
larger bosonic hard core diameter migrates to the edges 
of the trap. 

Shi et al. [l3| studied the phase separation of two- 
species trapped and untrapped Bose gases at finite tem- 
perature and found that the interspecies interactions af- 
fect the formation and depletion of the two condensates 
and lead to spatial phase-separation of the mixture. They 
argued that the shell is trapped in an effective potential 
which has a minimum away from the center of the trap 
close to the surface of the core. This effective trap is a 
combination of the traps confining the mixture and the 
interspecies interactions. According to Shi et al. then, 
condensation of the shell happens at the surface of the 
core and indeed we do observe a condensate in the shell 
as displayed in Figs. [Hand Figs. [51 
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C. Why do we use a HCSW? 

Particularly the HCSW is a suitable potential to de- 
scribe the attractive interactions between the HS bosons 
in this work, since it is a HC contact interaction plus 
an attractive tail added to the HC and the HC diameter 
is the same as the HS diameter of the bosons. Another 
reason for choosing the HCSW is to simulate a Feshbach 
resonance. This is because the HCSW has a well defined 
range, width, and depth and via these parameters one 
can easily tune the scattering length to be at the Fes- 
hbach resonance using Eq. ljlip below when a — s- ±00 
in order to check any instabilities (or stabilities) arising 
from this. In terms of designing the trial wave function, 
the exact solution of the two-body Schrodinger equation 
interacting via a HCSW lead us in the construction of a 
flexible Jastrow function for HS bosons with attractive 
interactions as mentioned in Sec. IH Dl 

D. HCSW parameters used 

Our main purpose for the choices of the previous val- 
ues of the HCSW parameters in this work was to provide 
a qualitative study of the properties of trapped Bose-gas 
mixtures with attractive interactions and to reach a qual- 
itative understanding of the role of the interatomic inter- 
actions in these properties. We first remind the reader 
that the values of the HC diameters have been chosen to 
enable substantial depletion of the condensate. 

1. Range 

The range of the HCSW {R = 0.59a;,o) used in this 
work is of the same order of magnitude as that used by 
Astrakharchik et al. [1^ for another model potential of 
the form V{r) = — Vb/ cosh^(r/ro). Here ro determines 
the range and they set vq = 0.1 Op where Cp is the trans- 
verse oscillatory trap length for a highly elongated trap. 
Their Op is small because of tight confinement along the 
transverse direction, similarly our a^o is also considered 
to be small since we use a tight trap. 

2. Depth 

We use a shallow HCSW which is much weaker than a 
realistic interatomic potential [3^. |40|. drawing our justi- 
fication from what has been noted before by Gao [4l| . He 
discussed improved interatomic model interactions be- 
yond the HS potential or delta function pseudopotential 
used in Gross-Pitaevskii theory. These model potentials 
are simple in a sense that they are shallow and are ap- 
plicable in quantum few-body and many-body systems. 

One of the most important points relevant to our work 
mentioned by Gao is that a real interatomic potential 
can simply become unmanagable if used in few-body or 



many-body calculations. A key conclusion in his paper 
is that the real potential in a many-body system around 
the threshold, such as a BEC state, no matter how deep 
this potential might be, can be replaced by an effective 
shallow potential that supports only one or two bound 
states. Gao shows that by using shallow model poten- 
tials, much weaker than the real potentials, the results 
are in good agreement with those using a real interatomic 
potential. 



E. Artificial stability of the mixtures, Feshbach 
resonance, and negative energies 

The HC potentials used in this investigation, whether 
attractive or repulsive, prevent real collapse and there- 
fore the mixtures are always stable and can not collapse 
to a singularity. In that sense, we speak of an "artifi- 
cial stability" This HC is present in the HCSW 
potential since upon "switching on" a HCSW one effec- 
tively adds to the HC potential an attractive tail, well 
defined in width and depth. Thus whatever the depth of 
the HCSW is, the bosons will not be able to approach 
each other to distances lesser than the HC diameter of 
the interactions (oc, be, or (ac + bc)/2) as imposed by the 
Jastrow functions. As mentioned in Section lllDi the Jas- 
trow function of the HCSW has a short-range repulsive 
and a long-range attractive part. The short-range part of 
the HCSW Jastrow keeps the bosons at some average dis- 
tance away from each other, whereas the attractive part 
tries to bring them closer together. The balance between 
the repulsive and attractive parts, keeps the system in 
equilibrium. In the case of repulsive interactions only, 
the system is primarily balanced by the repulsive HC 
and the external confining potential. 



1. At the Feshbach resonanee 

If one should increase the HCSW depths to values up 
to the first Feshbach resonance and beyond, the systems 
begin to shrink to very high densities. At this stage their 
energy is mainly potential (negative) and a large fraction 
of the bosons reside inside the HCSW. However they still 
show the artificial stability discussed above. For exam- 
ple in Fig[TH]we demonstrate how the stability of a two- 
species Bose gas of A^i = 20, = 10 particles, d = 0.54, 
and HS diameters Oc = 0.01, he = 0.02, respectively, 
shifts to higher a values as V12 is increased from shallow 
to deep values, even up to the first Feshbach resonance at 
Vi2 — —8.462 and beyond. The intraspecies interactions 
arc repulsive HCs. The HCSW depth corresponding to a 
Feshbach resonance is obtained from the condition that 
a ^ ±00, where a is the s-wave scatterin g len gth of the 
HCSW. And according to Giorgini et al. [43| a is given 

by 
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FIG. 18: Energy vs a showing stabilities at high HCSW depths 
for a two-component Bose gas of A'^i = 20 and N2 = 10 particles, 
ac = 0.01, be = 0.02, d = 0.54 and various interspecies HCSW 
depths V12 shown. The first Feshbach resonance for the system is 
at V12 = -8.462 



that. There is a cost in kinetic energy increase if the 
bosons move closer together. Thus the density of the 
hquid saturates eventually to a value set by the range of 
the HCSW (and less so by the HC diameter). 

There is a large energy release on condensation of the 
Bose gas to the dense liquid state. In FigllS] we show 
how the energy drops from E/N w 1.6 at V12 = to 
E/N « -90 at V12 = -10.0. Thus energy is released 
at a degree that is proportional to the well depth V. 
Upon condensation, there is substantial depletion of the 
condensate so that a large percentage of the atoms lie 
in states above the condensate. Thus the condensation 
process is characterized by large increase in the density of 
the system, a large drop in the energy of the system until 
there is a large release of condensation, and a substantial 
depletion of the condensate. 



F. Condensate depletion 



a = Rc + {R-Rc 



_ taii[KoiR - R,)] \ 



where Rc is the HC diameter (in our case Oc for com- 
ponent 1. be for component 2, or {ac + bc)/2 for the 
mixture), R the edge, and Kq = y/Vom/h? the wave 
vector of the HCSW. In trap units Rc Rc/aho = Rc 
and similarly for i?, Vq — > VQ/hioiio = Vq, and thus 

/Vo/aho 



Ko yVohhJhom/h? ^ VVo/aho = Kq. 

In Fig llSl at V12 = —10 the wave function is very much 
contracted and its density has risen substantially as in- 
dicated by a large Gaussian variational parameter 10. 
Note that one still obtains a deep negative energy min- 
imum at the first Feshbach resonance manifesting the 
strong stabilizing factor of the HCs. Thus even if we 
reach the Feshbach resonance and surpass it while vary- 
ing the HCSW depth, the systems remain artificially sta- 
ble. Going back to the previous Fig.|3] for example, we 
crossed a Feshbach resonance while increasing V, but no 
sharp density profile indicative of a collapse can be seen. 
In the next section, we shall explain the occurance of the 
negative energies to be the result of a liquefaction process 
of the Bose gas since when the energy becomes negative, 
energy is released from the system. 



2. Liquefaction of Bose gases (negative energies) 

The bosons can condenste to a dense liquid at critical 
values of the scattering length and potential depth. This 
is signified by the occurence of negative energies as in 
Figs.[2l[3l andflSl The liquid density is set by the range of 
the attractive well since there is a large energy incentive 
for a boson to lie in the attractive well of its neighbours 
but no further potential energy incentive to lie closer than 



In this section we explore possible reasons for the en- 
hanced depletion of the condensates in the mixtures. We 
believe that due to mixing the reduction in free volume 
between the HS bosons available for condensate forma- 
tion in each species is a common ground for enhanced 
depletion in all types of mixtures. The magnitude of re- 
duction in free volume varies however with the types and 
strengths of interactions. As the free volume between 
the HS bosons decreases, the probability for relocating a 
boson at a certain site, say ri and energy eo to another 
location at ri' and the same energy eo is reduced. This is 
because the chance of finding a site between the bosons 
large enough to accommodate a boson becomes lower. 

Going back to Fig. [7] then, the reason the mixtures 
with interspecies attraction (of Fig.|4]) indicated by (in- 
ter.) show a larger depletion for component 1 (with HC 
ac) than the mixtures with interspecies repulsion (of Fig. 
[5]) indicated by (intra.) is because the former are com- 
pletely mixed as compared to partial mixing of the lat- 
ter. In the case of complete mixing the available vol- 
ume for condensate formation is severely reduced and 
smaller than the case of partial mixing. Further, since 
both components in Fig. |4] are localized at the center of 
the trap, they contribute to their mutual condensate de- 
pletion where the density is highest, namely at the center 
of the trap. In the case of partial mixing the shell con- 
tributes to the depletion chicfiy at the edges of the trap 
and does not influence the condensate at the center of 
the trap very much. There could also be other reasons 
that explain the enhancement in the depletion. 

The scenario is however different for component 2. In 
Fig. [11] the depletion is larger for the 2BEC with inter- 
species repulsion (of Fig.[5|) than the 2BEC with inter- 
species attraction (of Fig. 2]). The reason may be be- 
cause the boson sizes Uc of the HC 2BEC are increased 
far above those of the HCSW 2BEC which remain fixed 
thus outweighing the role of the HCSW in the depletion 
(free- volume reduction) . 
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G. Ground state solutions 

Trippenbach et al. [l^ identified all possible classes 
of solutions for 2BECs and found that, in the case of 
isotropic harmonic trapping potentials, many spherically- 
symmetric phase-separated geometries are possible. In 
addition, symmetry breaking solutions do exist but 
within the TF approximation the ground state cannot 
be one with broken symmetry. Similarly, our mixtures 
are spherically symmetric in their ground states. 

V. CONCLUSIONS 

In summary then, we have investigated the effect of 
intra- and interspecies interactions on the properties of 
ultracold 2BECs in tight harmonic traps using VMC. The 
repulsive inter- or intraspecies interactions were mod- 
elled by a HC contact potential, the radius of which 
being equivalent to the s-wave scattering length in the 
low-energy and long-wavelength approximation. The at- 
tractive interactions were modelled by a HC repulsive 
part plus a shallow attractive well, the HCSW. We did 
not describe the attractive interactions by the HCSW 
scattering length a, but rather the depth of the HCSW 
in order to avoid the large fluctuations in the value of 
a. We calculated the energies, density distributions and 
condensate density distributions. We further obtained 
the condensate fractions of the components from the OB- 
DMs. A key point is that we chiefly focused on the role 
of interactions in enhancing the condensate depletion of 
each component in a mixture as compared to the case 
when each component is in a separate trap of its own. 
To the best of our knowledge, this has not been done in 
the previous literature on mixed Bose gases up to this 
date. We present novel physics associated with the effect 
of intra-and interspecies interactions on the condensate 
properties of mixed Bose gases. Wc find that: 

a) the mixing of two Bose gases in a trap enhances 
the condensate depletion of each gas as compared 
to the case when either one is in a separate trap of 
its own. In both cases of attractive and repulsive 
interactions the reduction in the available volume 
for condensate formation due to mixing plays a key 
role in the enhancement of the depletion. In the 
case of attractive interspecies interaction the en- 
hanced depletion may be further driven by the liq- 
uefaction of the Bose gases at the higher densities 
and the release of energy. 

b) when the condensates are phase separated due to 
strong repulsive interspecies interactions, the core 
remains stable and is not "squeezed" substantially 
by the shell. Complete mixing is still possible up 
to some repulsion threshold. 

c) according to Refs. [lB|, our mixtures are stable 
because they are spherically symixietric. 



d) we anticipate that the healing length of a IBEC 
changes upon mixing it with a second component 
into the system. 

c) in the case of complete phase-separation, although 
the density of the shell is much smaller than the 
core, substantial depletion is still observed in the 
shell triggered by the presence of large bosons in 
the shell. 

f) in the case of intercomponent attraction a 2BEC 
behaves similarly to a IBEC as the two compo- 
nents are completely mixed and allow the conden- 
sates of either component to migrate towards the 
edges of the trap at the higher densities. This is 
contrary to the somewhat opposite case of repul- 
sive interspecies interactions where the condensate 
of the core is pushed back towards the center of the 
trap. 

g) finally the HC potentials provide a strong stabiliz- 
ing mechanism for the Bose gases with attractive 
interactions. 
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APPENDIX A: MEAN-FIELD MODEL FOR THE 
ESTIMATION OF THE MIXTURE-ENERGIES 

We consider two mixed Bose gases of A^i and N2 par- 
ticles, HC diameters Oc and 6c, and bosonic masses mi 
and TO2, respectively, where initially the interspecies in- 
teractions are set to zero. That means the two Bose 
gases are initially independent of each other and both of 
them are concentric spheres at the trap center. We then 
construct a rough model that describes the energy of a 
boson-boson mixture by using the following assumptions. 
In our estimate for the energies, we derive our concepts 
from a paper by Ao and Chui [lH| who gave a simpli- 
fied expression for the total energy of an inhomogcneous 
binary Bose gas, 

(Al) 

where Gy = AttIt? aij / rriij are the interaction parame- 
ters, ai2 the inter- and an the intraspecies s-wave scat- 
tering lengths, ma the mass of a boson in one compo- 
nent and mi2 the reduced mass and V is the volume 
of the gas. This expression neglects the kinetic energy 
(quantum pressure) of each component. We modify this 
expression by replacing the first two terms on the right- 
hand-side by the Thomas-Fermi (TF) energy of each com- 
ponent. We use trap units in terms of component 1, i.e 



19 



aho = \/ h/miLUi and hoji for both systems as done be- 
fore. The TF energy for each component is then 



(A2) 



where i = 1, 2, /^i ~ i (15 A^i Uc)^^^ is the chemical po- 
tential of component 1 and 



fJ-2 



2 



1/5 



(A3) 



that of component 2. The TF radius of component 1 is 
Rtfa = (ISA^iCc)^/^ and that of component 2 



Rtf.2 = [IbN^Kf/" 



2/5 



(A4) 



Imagine now switching the interspecies interactions on 
such that the Bose gas with the larger HC diameter is 
expelled towards the edges of the trap and forms a shell. 
The shell would then lie approximately at the TF radius 
of the core away from the center of the trap. The vol- 
ume of the cloud is then approximately y = 4 tt R^p ^/3 

if Oc < bc&wd AttR^p^/-^ if "^c > ^c- Thus the size of 
the cloud in largely determined by the core in the case 
of moderate repulsive interactions. In order to calcu- 
late its energy, we therefore assume a superposition of 
its initial TF energy when both mutually noninteracting 
components are localized at the center of the trap and 
an approximate potential energy for the particles of the 
shell formed at the edges of the trap after switching on 
the interspecies interactions. As a result, component 2 
gains additional potential energy beyond Etf,2 when it 
becomes a shell. An estimate for the potential energy of 
the shell is 



Vtrap = ^NimiUj1^R^p2/fnOho 

4/5 



\n, ( ^ ) {l^N.h^f" (A5) 

2 V"l2 / 



if component 1 forms a shell and 



Vtrap = -^N2m2Ujl„RTp^i/fuJho = 

In, f (15N,a^f' (A6) 
2 \mi J 



if component 2 forms a shell. That is we used 



R) 



'TF,2 



and 



EN2 
^=l ^2 - 

«=1 ™2 



2i 



{r,} ~ Rtfa 



(A7) 



(A8) 



in estimating the radius of gyration of the core. Note 
that mk ~ Nknik {k — 1 or 2) is the total mass 

of either component as all the particles in a component 
have the same mass. The interspecies interactions can be 
calculated according to lAll 



V GiiG 
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V 



(A9) 



Gathering all the previous terms together, the total en- 
ergy of the mixture is then 



E 



MFA 



{Ni + N2) 
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iVi(15iVia,)2/5 + iV2(15iV2&c)'/'(TO2/™l)'/' 



■TF,1 
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nrppo 



R 



0-c < be 
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iiVi (toi/to2)^/' (15iV2&c)"/' : a,>b, 
iiV2 (m2/mi) (15iViaj2/5 ; a, < b^ 



(AlO) 



where as a reminder oji =102- 
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